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We report a new approach to simulate amorphous networks of covalently bonded materials that 
leads to excellent radial distribution functions and realistic atomic arrangements. We apply it to 
generate the first ab initio structures of nitrogen-doped silicon, a-SiN^;, for thirteen values of x from 
to the nearly stoichiometric composition of x=1.29, using the Harris functional and thermally 
amorphizised, periodically continued, diamond-like cells with 64 atoms. Partial radial features are 
reported for the first time and the total radial distribution functions agree very well with the few 
existing experiments. Our results should stimulate further experimental and theoretical studies in 
amorphous covalent materials. 

PACS numbers:71.23.Cq, 71.15.Pd, 71.55.Jv, 

I. INTRODUCTION 



The properties of amorphous silicon-nitrogen alloys 
have attracted a great deal of attention in the last decade; 
a-SiNj; has electrical, optical and mechanical features 
useful in a variety of industrial applications and the 
strong covalency of their bonding makes them the pro- 
totype of covalent materials. Their optical gaps depend 
strongly on the nitrogen content a; for < x < 1.33 so 
they can be tuned to fit specific needs in solar cells; their 
total and partial radial distribution functions (RDFs) are 
practically unknown, except for the stoichiometric con- 
tent. Some semiempirical studies have been done on their 
electronic structure, optical gaps and RDFs and a first- 
principles approach has been used on an ad /loc-generated 
amorphous structure. Therefore, any ab initio approach 
that adequately generates, describes and predicts fea- 
tures of a-SiN^: may have a wider applicability to deal 
with other covalently-bonded amorphous solids; in par- 
ticular, systems like a-GeN^; should be amenable to our 
approach [Q. 

Recently [| we carried out ab initio Harris-functional- 
based studies of the atomic and electronic structure of 
pure and hydrogenated amorphous silicon, using 64 sil- 
icon atom cells plus hydrogens that diffuse within the 
cells. We now apply these methods to amorphous silicon- 
nitrogen alloys to test the adequacy and the predictive 
powers of our approach for covalent materials. To the 
best of our knowledge, these are the first ab initio ther- 
mally generated amorphous networks where studies of 
the atomic topology of a-SiN^ are carried out; the con- 
tent X is in the range < a; < 1.29 where x — y/{Q4: — y) 
and y is the number of nitrogen atoms. 

The experimental and theoretical activity prior to 1990 
is well documented in a paper by Robertson where per- 
tinent references can be found In this work Robert- 
son reports his semiempirical tight-binding calculations 
for the optical gaps and also reports several experimen- 
tal gaps for hydrogenated and non-hydrogenated o-SiN^; 
alloys. Bethe lattice calculations have been done by 



Marti'n-Moreno et al. ^ and by San-Fabian et al. ||5| also 
using semiempirical parameters, whereas Ordejon and 
Yndurain |^ in a very nice work do non-parameterized 
calculations of a-SiN^; where the equilibrium positions of 
Si and N atoms in clusters are ported to the alloy network 
constructed ad hoc. They obtain a wealth of information 
including optical gaps, however tetrahedral coordination 
of the silicon atoms and threefold planar coordination 
of the nitrogen atoms is assumed with interatomic dis- 
tances of 2.33 A for Si-Si and 1.74-1.76 A for Si-N. A 
general characteristic of these calculations/simulations is 
that gap states, when considered, are introduced either 
by hand, progressively replacing Si by N, or by the al- 
gorithms that generate the random networks, unlike the 
procedure reported here. Recent semiempirical classical 
simulations by de Brito et al. [Q produced total RDFs 
and average nearest-neighbor results that are the subject 
of comparison with ours. 

On the experimental side, as early as 1976 
Voskoboynikov et al. ^ studied some RDFs and optical 
gaps of silicon-rich silicon-nitrogen films as a function of 
the gas ratio. It was then observed that the gaps increase 
as a function of the nitrogen content; the films seemed to 
contain hydrogen and large clusters of silicon. However, 
reliable experimental RDFs are scarce Q and, appar- 
ently, only total ones for the stoichiometric amorphous 
composition exist |l^ , where a decomposition of the sec- 
ond peak of the total RDF into its partial contributions 
is also carried out. 

It is clear that the atomic topology also determines 
the electronic properties of the amorphous samples, and 
therefore any understanding of the RDFs and the atomic 
distribution in the random networks is relevant in the 
characterization of the electronic and optical properties 
of these materials. In what follows we report the genera- 
tion of random networks for amorphous silicon nitrogen 
alloys that leads to RDFs in good agreement with what 
is known experimentally and predict atomic structures 
for a variety of nitrogen contents. Electronic gap states 
are found in these structures that shall be dealt with in 
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a future publication. 



II. METHOD 

Our 26 amorphous samples of a-SiN^; were generated 
with FastStructure a DFT code based on the Har- 
ris functional produced by Molecular Simulations, Inc.. 
The optimization techniques use a fast force generator 
to allow simulated annealing/molecular dynamics stud- 
ies with quantum force calculations |12|. The LDA pa- 
rameterization invoked is that due to Vosko, Wilk and 
Nusair ||l3|. An all electron calculation is carried out, 
and a minimal basis set of atomic orbitals was chosen 
with a cutoff radius of 5 A for the amorphization and 
3 A for the optimization. The physical masses of nitro- 
gen and silicon are used throughout and this allows us to 
see realistic randomizing processes of the atoms during 
the amorphization of the supercell. Finally, the forces 
are calculated using rigorous formal derivatives of the 
expression for the energy in the Harris functional JT^ . 

In order to test the adequacy of calculations carried 
out with FastStructure we used it to obtain the size of 
the crystalline cell of /3-Si3N4 that minimizes the energy 
at the F-point. Fig. 1 shows the results of such cal- 
culation. The experimental crystalline volumen is given 
by 145.920 A'^ whereas the calculated volumen is 
146.797 A"^; a deviation of 0.6%. For this reason we feel 
cautiously optimistic about the use of FastStructure to 
generate random networks of silicon-nitrogen alloys. 

Since it has become increasingly clear that quenching 
from a melt generates undesirable structures Q we took 
a different path We amorphisized the crystalline di- 
amond structures with a total of 64 atoms ((64 — y) sili- 
cons and y substitutional nitrogens) in the cell by slowly 
heating it, linearly, from room temperature to just be- 
low the corresponding melting point for each x, and then 
slowly cooling it to K. To determine the melting tem- 
peratures for each x we linearly interpolated between 
the pure silicon value and the stoichiometric compound, 
X = 4/3 = 1.33, and then remained below these temper- 
atures (Table I). Since the time step was the same for all 
runs, 6 fs, and the melting temperatures increased with 
x, the heating/cooling rate varied from 2.30 x lO^^K/s for 
pure silicon, to 3.11 x lO^^K/s for x = 1.29. The atoms 
were allowed to move freely within each cell whose vol- 
ume was determined by the corresponding density and 
content. Table L The densities were taken from the ex- 
perimental results of Guraya et al. Once this first 
stage was completed, we subjected each cell to annealing 
cycles with intermediate quenching processes. Finally 
the samples were energy-optimized to make sure the fi- 
nal structures would be those of a local energy minimum. 



TABLE I 

Contents, Melting Temperatures and Densities 
for a-SiNr 



Sample 


X 


Melting 


Density 






Temp. (K) 


(g/cc) 




0.000 


1680 


2.329 




0.085 


1747 


2.435 


Si Ni n 


0.185 


1814 


2.512 


Si/ioNi c; 


0.306 


1881 


2.600 


l)1 a a in on 


0.455 


1948 


2.694 


Si39N25 


0.641 


2015 


2.803 


Si34N30 


0.882 


2082 


2.931 


Si33N31 


0.939 


2095 


2.957 


Si32N32 


1.000 


2108 


2.988 


Si3lN33 


1.065 


2122 


3.017 


Si3oN34 


1.133 


2136 


3.048 


Si29N35 


1.207 


2149 


3.081 


8128^36 


1.286 


2162 


3.115 



It should be kept in mind that our objective is always 
to generate realistic amorphous structures of a-SiNj. and 
not, in any way, to mimic the experimental processes 
used to produce these alloys. 



III. RESULTS AND DISCUSSION 

We performed two runs for each x value and from 
X = 0.882 on (a-Si34N3o) the number of nitrogens was 
increased one at a time to be able to map the interesting 
processes that occur for these contents (percolation of 
the Si-Si bonds, widening of the optical gaps, etc). Once 
the atomic structures were obtained, we calculated the 
corresponding total and partial RDFs for each of the 26 
runs and averaged them by corresponding pairs. Of those 
13 averaged plots we are reporting here total and par- 
tial RDFs for pure silicon, for the almost stoichiometric 
sample, a-Si28N3g, x = 1.29, and two intermediate ones: 
a-Si44N2o, X = 0.46 and a-Si32N32, x = 1.00. Fig 2a for 
the pure amorphous silicon sample shows our results and 
the upper and lower bounds of the available experimen- 
tal data the agreement is good since our RDF falls 
within these bounds and the four experimental peaks are 
correctly reproduced by our simulations. Figs. 2b to 2d 
show the variation of the partial RDFs for Si-Si, Si-N 
and N-N as a function of content and their contribution 
to the total RDF. In Fig. 2d the composition of the sec- 
ond peak of the total RDF for the nearly a-Si3N4 can be 
observed and it agrees completely with experiment 
since it is formed by the average second-neighbor (2n) 
contributions of mainly the N-N and Si-Si partials and 
to a lesser extent by the Si-N partial. The third peak is 
essentially due to the Si-N partial with a small contribu- 
tion from the N-N partial but, as far as we know, there 
are no experimental results for comparison. 
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Figs. 2 show that as the nitrogen content increases 
the first peak of the total RDF (f .85 A), which is due to 
the Si-N average nearest-neighbor {(nn)) contributions, 
increases systematically and the (nn) Si-Si peak (2.45 A) 
decreases systematically. The third peak moves toward 
low r as X increases since the N-N contribution becomes 
more predominant at high content (3.25 A to 2.95 A), fn 
our structures there are no (nn) nitrogens since the con- 
tent is below stoichiometry and nitrogens have a marked 
tendency to bind to silicons. For the nearly stoichiomet- 
ric sample, x = f .29, the Si-Si (nn) contribution to the 
total RDF has practically dissapeared and this implies 
that there is a nitrogen atom between every pair of sili- 
cons. This is borne out by the results presented in Figs. 

3 where a study of the the average coordination numbers 
(cn) in the 13 random networks is depicted. The follow- 
ing cutoff radii were used: Si-Si, 2.55 A; N-N, 3.35 A; 
and Si-N, 2.15 A, which are the positions of the minima 
after the first peaks of the corresponding partials. Fig. 
3a shows the results of our simulations where it can be 
seen that the N-N plot flattens for x « 1.1, the percola- 
tion threshold of Si-Si bonds ||l^ ; the Si-Si (nn) go from 

4 to practically 0. The Si-N graph refers to the (nn) ni- 
trogens around the silicon atoms and varies from to 4, 
whereas the N-Si refers to the (nn) silicons around nitro- 
gens and indicates that nitrogens immediately surround 
themselves with practically 3 Si, saturating its valence. 
The crossing of the Si-Si and Si-N plots at a; « 0.7, is in 
agreement with experiment (Davis et al. |^). There is a 
crossing of the Si-Si, N-Si and N-N plots at x w 0.3 and 
a crossing of N-Si and Si-N at a; w 1.0 which have been 
observed experimentally for hydrogenated alloys by Gu- 
raya et al. Fig. 3b. However, due to the presence of 
hydrogen a curvature appears in the (cn) for Si-Si, Si-N 
and N-Si so, in order to compare our results to this ex- 
periment, we did the following. We carried out the sum 
of N-H plus N-Si from the experiment, the average total 
number of atoms that surround a N, N-*, and plotted it 
along with our N-Si; we also did the sum of the experi- 
mental Si-N, Si-Si and Si-H, the average total number of 
atoms that surround a Si, Si-*, and plotted that along 
with our sum of Si-N plus Si-Si. This is presented in Fig. 
3c. It is clear that our predictions closely agree with 
the integrated experimental results and validate our ap- 
proach. The discrepancies are most likely due to existing 
dangling and floating bonds. 

Fig. 4 is the comparison of our ah initio results with 
the classical Monte Carlo simulations of de Brito et al. 

where they used empirical potentials developed a la 
Tersoff for the interactions between Si and N. It is clear 
that although the positions of some peaks are reproduced 
in both simulations, the general behavior of the total 
RDFs only agree qualitatively. Comparison of each of 
these simulations with experiment is presented in Fig. 5 
where agreements and discrepancies can be appreciated. 



IV. CONCLUSIONS 

We have devised an ah initio approach for a-SiN^, 
(0 < X < 1.29) that generates radial distribution func- 
tions in agreement with existing experimental results and 
realistic atomic structures where the simulated average 
coordination numbers and the experimental ones coin- 
cide. Total RDFs agree very well with experiment, where 
available, and partial RDFs show that the Si-Si (nn) peak 
disappears as nitrogen increases indicating a tendency to 
form 6-atom arrangements as the content x approaches 
the stoichiometric value. Experiment shows that for a- 
Si3N4 Si and N form closed rings, Si-N-Si-N-Si-N, typical 
of the crystalline Si^Ni structures. The growth of the 
Si-N peak as nitrogen increases bears out this behavior. 
No nitrogen- nitrogen bonds, including molecular nitro- 
gen, are observed in the final structures even though for 
X > 1 the starting diamond structure does contain nitro- 
gens next to one another. For x ^ 1.1 the effects of the 
percolation threshold of the Si-Si bonds is observed in 
the N-N (2n). For x « 0.7 the Si-Si and Si-N neighbors 
are practically the same, as found experimentally. Also, 
Si-Si, N-Si an N-N are practically the same for x w 0.3 as 
are Si-N and N-Si for x w 1.0. The integrated experimen- 
tal results and our simulation agree. The first prominent 
peak in the total RDF of the nearly stoichiometric sam- 
ple is due to Si-N and an analysis of the second peak 
indicates that N-N, Si-Si and Si-N contribute to it, in 
agreement with experiment. The third peak is mainly 
due to the Si-N, with a small contribution from the N- 
N; no experimental results exist for comparison. Our 
approach, being ah initio, is of wider applicability than 
classical or semiempirical ones and should be relevant for 
the understanding of the physics of amorphous covalent 
materials. 
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Figurel. Energy vs volume for crystalline beta-Si3N4 
obtained with FastStructure and experimental values 
taken from Wild et aL (see text). The experimental 
and calculated volumes agree to within 0.6%. 

Figure2. Total and partial RDFs for (a) pure silicon; 
the light curves are the upper and lower experimental 
bounds (Rcf. [2]); (b) a-Si44N20, x=0.46; (c) a-Si32N32, 
x=1.00; (d) the nearly stoichiometric sample a-Si28N36, 
x=1.29. 

Figure3. Average coordination numbers jcn/^ as a func- 
tion of X. (a) Our results, (b) Experimental results for 
hydrogenated alloys from Guraya et al.. (c) Comparison 
of the integrated results (see text). Lines are drawn as 
guides to the eye. 

Figure4. Comparison of our results and those obtained 
by de Brito et al. (Rcf. [7]) for the total RDF. The 
agreement is at best qualitative although the position of 
some of the peaks coincide. 

Figure5. (a) Comparison of our simulations and ex- 
periment (Ref. [10]) for the stoichiometric sample, (b) 
Comparison of de Brito's simulations and experiment for 
the same sample. 
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